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ABSTRACT 

Free  surface  flows  of  a  liquid  poured  from  a  container  are  calculated 
numerically  for  various  configurations  of  the  lip.  The  flow  is  assumed  to  be 
steady,  two  dimensional  and  irrotational,  the  liquid  is  treated  as  inviscid 
and  incompressible,  and  gravity  is  taken  into  account.  It  is  shown  that  there 
are  flows  which  follow  along  the  under  side  of  the  lip  or  spout,  as  in  the 
well-known  "teapot  effect*,  which  was  treated  previously  without  including 
gravity.  Some  of  the  results  are  applicable  also  to  flows  over  weirs,  and 
spillways . 
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SIGNIFICANCE  AND  EXPLANATION 


Wie  flow  of  a  liquid  poured  from  a  container  is  a  free  boundary  flow 
driven  by  gravity,  and  it  has  proved  difficult  to  determine  such  flows 
analytically.  Therefore  we  have  developed  a  method  to  calculate  them 
numerically. 

Two  calculated  pouring  flows  from  a  thin-walled  vessel  are  shown  in 
Figures  1  and  2.  The  flow  in  Figure  1  is  that  desired  when  pouring  a  beverage 
like  tea,  while  the  flow  in  Figure  2  is  the  undesirable  one  which  occurs  in 
the  "teapot  effect". 

Flows  which  run  along  a  wall  are  desirable  when  pouring  beer  to  avoid 
excess  foam,  and  when  pouring  acid  to  avoid  splashing.  All  these  flows  are 
like  the  calculated  flow  shown  in  Figure  3,  which  also  represents  the  flow  in 
a  spillway. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


POURING  FLOWS 


Jean-Marc  Vanda n-Broeck ^  and  Joseph  B.  Keller  ' ^ 


1.  Introduction 

The  flow  of  a  liquid  poured  from  a  container  is  a  free  boundary  flow  driven  by  gravity,  and  it  has 
proved  difficult  to  determine  such  flows  analytically.  Therefore  we  have  developed  a  method  to 
calculate  them  numerically  when  the  flow  is  steady,  two  dimensional  and  irrotational,  and  the  fluid 
is  inviscid  and  incompressible.  We  have  used  it  before  to  calculate  flows  over  weirs  [1],  and  now  we 
shall  adapt  it  to  the  calculation  of  flows  over  lips  or  spouts  of  various  shapes.  Some  of  our  results 
are  applicable  also  to  weirs  and  spillways. 

Two  calculated  pouring  flows  from  a  thin-walled  vessel  are  shown  in  Figures  1  and  2.  In  both 
cases  the  angle  0  between  the  vessel  wall  and  the  horizontal  is  0  =  x/3.  The  flow  in  Figure  1  is 
that  desired  when  pouring  a  beverage  like  tea,  while  the  flow  in  Figure  2  is  the  undesirable  one 
which  occurs  in  the  uteapot  effect”.  It  was  analyzed  previously  without  taking  account  of  gravity 
[2],  and  the  present  calculation  shows  that  such  a  flow  can  occur  also  when  gravity  is  taken  into 


Figure  2.  The  undesirable  "teapot  effect"  flow  over  the  same  thin  wall 

or  weir  as  in  Figure  1.  The  flow  shown  has  one  free  streamline 
It  was  calculated  for  8  =  tt/3  by  the  second  method  in  Section 
5.  The  x,y  coordinates  are  shown  with  their  origin  at  C, 
the  end  of  the  wall.  The  points  I  and  J  are  at  infinity 


on  opposite  sides  of  the  wall. 


account. 


A  flow  like  that  in  Figure  1,  but  with  $  =  x/2,  was  calculated  in  [1]  to  describe  the  flow  over 
a  thin  vertical  weir  in  a  deep  channel.  In  the  same  way,  the  flow  in  Figure  1  applies  to  a  thin  weir 
sloping  at  the  angle  /3  —  t/3  in  a  deep  channel.  Flows  for  other  slopes  can  be  calculated  by  the 
present  method. 

Flows  which  run  along  a  wall  are  desirable  when  pouring  beer  to  avoid  excess  foam,  and  when 
pouring  acid  to  avoid  splashing.  But  then  the  flow  must  run  along  the  inner  wall  of  the  receiving 
container,  rather  than  along  the  outer  wall  of  the  pouring  container.  In  chemical  laboratories  to 
avoid  splashing,  liquids  are  often  poured  onto  stirring  rods  which  they  run  along.  All  these  flows 
are  like  the  calculated  flow  shown  in  Figure  3,  which  also  represents  the  flow  in  a  spillway.  For  it 
the  angle  between  the  two  walls  is  /?  =  3jt/4  and  the  Froude  number  F  =  U(gH)~l/2  is  F  =  1.3. 

Similar  flows  can  be  calculated  for  any  F  >  1  and  for  any  value  of  /?  satisfying  0  <  /?  <  x. 
Another  example,  shown  in  Figure  4  for  /?  =  x/4  and  F  =  1.3,  represents  the  undesirable  flow  over 
a  lip  which  has  its  sides  meeting  at  the  angle  /?.  Goh  and  Tuck  [3]  have  calculated  the  desirable 
flow  from  a  spout  consisting  of  two  horizontal  plates. 

In  Section  2  we  shall  formulate  the  problem  for  the  flow  shown  in  Figure  4  for  any  F  >  1 
and  /?  <  x ,  and  we  shall  present  the  analytical  solution  of  this  problem  for  F  =  oo  and  show  how 
to  obtain  the  asymptotic  form  of  the  solution  for  F  ^  1.  In  Section  4  we  describe  our  numerical 
results,  which  yield  the  free  streamlines  shown  in  Figures  3  and  4.  In  Section  5  we  indicate  how  to 
modify  the  method  presented  in  Section  2  in  order  to  treat  the  flows  shown  in  Figures  1  and  2. 

2.  Formulation 

Let  us  consider  the  flow  in  the  region  of  the  z- plane  shown  in  Figure  4.  The  horizontal  bottom  1C  is 
a  streamline  on  which  we  require  that  the  stream  function  V  =  0.  The  wall  CJ,  which  slopes  at  an 
angle  (i  from  the  horizontal,  is  part  of  the  same  streamline.  We  choose  cartesian  coordinates  with 
the  y-axis  along  the  horizontal  wall  IC  and  the  x-axis  directed  vertically  downwards.  The  origin  is 
at  the  corner  C  and  gravity  acts  in  the  positive  x  direction.  As  x  — +  +oo,  the  flow  approaches  the 
thin  jet  flow  of  Keller  and  Weitz  [4]  and  Keller  and  Geer  [5].  As  y  —*  —  oo,  the  flow  approaches  a 
uniform  stream  with  a  constant  velocity  U  in  the  y-direction.  The  free  surface  I J  is  a  streamline 
on  which  tp(x,y)  =  UH.  Here  H  is  the  depth  of  the  fluid  at  y  =  —  oo. 

Let  the  complex  potential  be  /  =  +  »0.  Without  loss  of  generality  we  choose  <p  =  0  at  the 

corner  C.  On  the  free  surface  IJ  where  the  pressure  is  constant,  the  Bernoulli  equation  yields 


the  asymptotic  solution 


I 

i/«* 


Figure  4.  The  "teapot  effect"  flow  over  a  wedge-shaped  lip  of  angle  6. 

The  points  I  and  J  are  at  infinity  on  opposite  sides  of 
the  wedge.  The  flow  along  the  horizontal  wall  has  thickness 
H  and  velocity  U  at  I.  The  x,y  coordinates  have  their 
origin  at  C,  the  vertex  of  the  wedge,  with  x  increasing  in 
the  downward  direction  and  y  increasing  horizontally  to  the 
right.  The  streamline  shown  was  computed  by  the  method  of 
Section  2  with  F  =  1.3  and  0  =  tt/4.  The  broken  line  is 


computed  from  the  asymptotic  solution  (3.6). 


''M 
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We  choose  H  as  the  unit  of  length  and  V  as  the  unit  of  velocity.  Then  (2.1)  becomes 


2_  i  2 


-  j^x  =  1  +  —  on  ip=  1. 


Here  F  is  the  Froude  number  defined  by 

F  =  V/(gH)''2.  (2.3) 

We  shall  restrict  our  attention  to  supercritical  flows,  for  which  F  >  1.  The  plane  of  the  complex 
potential  /  =  <p  +  iif>  is  shown  in  Figure  5a. 

Let  the  complex  velocity  be  {  =  u  —  tv,  where  tt  and  v  are  the  x  and  y  components  of  the 
velocity.  For  F_1  0  (i.e.,  g  ^  0),  the  velocity  £  increases  like  /*/3  as  ip  — *■  oo  [4].  Thus  we  have 

C  ~  /1/3  as  <p -*■  +oo,  F-1  ^  0.  (2.4) 

At  <fi  =  — oo  the  flow  is  supercritical  and  is  characterized  by  the  presence  of  exponentially  decreasing 
terms  in  £.  Thus  (  has  the  form 

^  —  *[1  4-  Ae*x/]  as  <p  — ►  -oo.  (2.5) 

Here  A  is  a  constant  to  be  found  as  part  of  the  solution  and  A  is  the  smallest  positive  root  of 

*■  A  —  F~2  tan  xA  =  0.  (2.6) 

Near  C,  where  there  is  a  corner  of  angle  /?,  the  flow  has  the  corner  behavior 

< A  as  /  — ►  0.  (2  7) 

The  problem  is  to  find  (  as  an  analytic  function  of  /  =  <p  +  iij>  in  the  strip  0  <  0  <  1,  satisfying 
(2.2),  (2.4),  (2.5),  (2.7)  and  the  kinematic  conditions 

u  =  0  on  V”  =  0,  <p  <  0,  (2.8) 


u  = -v  tan/9  on  V’  =  0,  ip  >  0. 


We  define  the  new  variable  t  by  the  relation 


/=*ln i±i 

J  X  1  -  t 


(2.10) 


v\  s 


The  transformation  (2.10)  maps  the  flow  domain  into  the  interior  of  the  onit  circle  in  the  t- plane 
so  that  the  walls  go  onto  the  real  diameter  and  the  free  surface  goes  onto  that  portion  of  the 
circumference  lying  in  the  upper  half  of  the  t  plane.  (See  Figure  5b). 

We  define  the  function  Q(t)  by  the  relation 

C  =  -i(-t)^/'_,[-lne(l  -<)]1/3[-ln2c]1/3[l  +  (1  +  t)2Aft(t)].  (2.11) 


Here  c  is  a  real  constant  between  0  and  1/2.  We  shall  choose  c  =  0.2.  It  can  be  checked  easily  that 
the  expression  (2.11)  satisfies  the  conditions  (2.4),  (2.5)  and  (2.7).  The  function  Q(t)  is  analytic 
for  |f|  <  1  and  continuous  for  jf|  <  1.  The  kinematic  conditions  (2.8)  and(2.9)  imply  that  the 
expansion  of  fl(t)  in  powers  of  t  has  real  coefficients. 

With  this  expansion  inserted,  (2.11)  becomes 


C  = 


-i(-*yv,r-1[-  In c(l  -  t))1/3[-  In 2c]-1/3 


l  +  (l  +  t)2A£antn 

n=0 


(2.12) 


For  given  values  of  F  and  /?,  the  unknown  coefficients  a„  have  to  be  determined  to  make  (2.12) 
satisfy  the  Bernoulli  condition  (2.2). 

We  use  the  notation  t  =  \t\eio,  so  that  points  on  the  free  surface  are  given  by  t  =  e‘er,  0  <  a  <  r. 
We  find  it  convenient  to  eliminate  *  from  (2.2)  by  differentiating  it  with  respect  to  a.  This  yields 
- =  0.  Now  on  the  circular  arc  f  =  ln(cotf).  Differentiating  this  relation  and  using 


the  identity 

dx  ^  dy  _  1 
d<p  *  dffi  C 


(2.13) 


we  obtain 

dx  _  _ 1 _ u(<y) 

da  ~  2  sin  a  [u(<r)]J  +  (w(<7)]2 

dy  1 _ v(a) 

da  ~  2  sin  a  (u(<r)]2  +  [w(<x)]2 


(2.14) 

(2.15) 


Here  ((a)  =  u(<r)  -  iv(a)  denotes  the  value  of  C  on  the  free  surface.  Upon  inserting  (2.14)  into  the 
differentiated  form  of  the  constant  pressure  condition,  we  obtain 


«(cr)ti«r(<r)  +  v(a)v„(a)  + 


r2  p(<r) 

2  sin  a  [ti(cr)]2  +  [r(<r)]2 


=  0 


(2.16) 


We  now  set  t  =  e'°  in  (2.12)  to  get  £(<r)  and  substitute  that  expression  into  (2.16).  We  will 
use  the  resulting  equation  to  get  the  unknown  coefficients  on.  To  do  so,  we  truncate  the  infinite 


series  in  (2.12)  after  N  terms.  We  find  the  N  coefficients  a.,  n  =  0, . . . ,  N  —  1  by  collocation.  Thus 
we  introduce  the  N  mesh  points 

+  /  =  1*-  (*•”) 

By  using  (2.12)  we  obtain  C(<r/)  and  d((<ri)/da  in  terms  of  the  coefficients  a„.  Substituting  these 
expressions  into  (2.16)  at  the  points  o/,  we  obtain  N  nonlinear  algebraic  equations  for  the  N 
unknown  a„.  This  system  is  solved  by  Newton’s  method. 

When  the  coefficients  an  are  known  for  a  given  F,  the  values  of  x(<r)  and  y(<r)  on  the  free 
surface  are  obtained  by  integrating  (2.14)  and  (2.15): 

i(<r)  =  *o  ~\J  S(<0[“atf(fi(fr)2  + 

y(<r)  =  fio  “  \  J*  e(<r)[6ino(u(«r)J  +  6(<r)2)]“1do 

Here  zo  end  yo  are  the  values  of  z  and  y  at  <p  =  0,  V>  =  1.  These  values  are  obtained  by  integrating 
dx/difr  and  dy/dij)  along  the  equipotential  <p  =  0  from  =  0  (i.e.,  from  the  origin)  to  =  1. 

S.  Free  streamline  solutions 

Before  describing  our  numerical  results  we  shall  show  that  the  problem  has  an  explicit  solution 
when  F  =  oo.  When  F  —  oo,  the  Bernoulli  condition  (2.2)  reduces  to 

ICI2  - 1  on  *=1.  (3.1) 

Then  (2.4)  must  be  replaced  by  the  requirement  that  £  remains  bounded  as  ip  — ►  +oo. 

It  can  be  checked  easily  that 

C  =  -i(-t)’-1  (3.2) 

satisfies  this  requirement  and  also  satisfies  the  conditions  (3.1),  (2.7)-(2.9).  Therefore  (3.2)  is  the 
solution  of  the  problem  for  F  =  oo.  Now  we  eliminate  f  between  (2.10)  and  (3.2)  and  use  the 
identity  (2.13).  Then  we  obtain  after  integration 

z  =  ijj  J-tanh^  **  df  (3.3) 

This  is  the  representation  of  the  solution  with  z  given  as  a  function  of  /.  For  0  =  0,  (3.3)  reduces 
to 

z  =  i  J  tanh  = -*~lncoth 


-11- 


A~  --  1 

*■ 


This  is  the  solution  derived  by  Keller  [2].  [See  his  equation  (10).] 

A  solution  for  F  large  can  be  obtained  by  considering  (3.2)  as  the  first  term  of  an  inner 
expansion  of  {(<,  F)  in  powers  of  F~2 .  This  expansion  is  to  be  matched  to  an  outer  expansion 
with(2.5)  as  leading  term  as  <p  — ►  — oo  and  to  another  with  leading  term  (2.4)  as  ip  — ►  +oo. 

By  nsing  (3.2)  and  (2.10)  we  find  that 

C  ~  -ijl  +  ^  -  l)c"^j  46  (3.4) 

Relation  (2.6)  shows  that  A  — >  1/2  as  F  — >  oo.  Therefore  matching  (2.5)  with  (3.4)  yields 

F  —  oo.  (3.5) 

As  <p  — ►  +oo,  the  solution  is  described  by  the  classical  thin  jet  theory  [5].  If  6  denotes  the  thickness 
of  the  jet,  we  have  |£|  ~  1/6.  Let  us  write  the  equation  of  the  free  surface  of  the  jet  as  y  =  y(x). 
Then  6  =  x  cos/3  +  y(z)  sin  0  ~  |C1— 1  -  Substituting  these  expressions  into  (2.2)  and  solving  for  y{x) 
yields 


y  =  y(x)  ~  -x  cot  0  +  - t—z{F 2  +  2  +  2x)~1/2. 

srn/3 


(3.6) 


Relation  (3.6)  is  valid  for  all  F  provided  that  x  >  F2.  An  easy  calculation  with  (2.2),  (2.13)  and 
(3.6)  shows  that  the  first  term  in  (3.6)  (i.e.,  y  ^  x  cot  0)  yields  (2.4). 

The  solution  (3.2)  gives  the  following  expression  for  the  shape  of  the  free  surface. 

1 


y  ~  -xco\0  + 


sin  0 


as  (fi  -*  -f  oo 


(3.7) 


Therefore  (3.7)  matches  with  the  leading  term  in  (3.6)  for  F  large. 

We  have  shown  that  (3.2)  matches  with  (2.4)  and  (2.5)  as  <p  — *•  +oo  and  — *  — oo  respectively. 
By  continuing  each  expansion  to  higher  order  in  F~2  and  matching,  wp  can  obtain  further  terms 
in  the  expansions.  Together  they  will  represent  the  flow  everywhere  for  F  large. 


4.  Numerical  results 


The  numerical  procedure  described  in  Section  2  was  used  to  compute  solutions  for  various  values  of 
F  and  0.  Solutions  for  0  <  r/2  illustrate  the  teapot  effect  whereas  solutions  for  0  >  x/2  represent 
flows  in  a  spillway.  Typical  profiles  with  F  =  1.3  for  0  -  3r/4  and  0  =  x/4  are  shown  in  Figures 
3  and  4  respectively.  The  broken  line  corresponds  to  the  asymptotic  solution  (3.5)  for  x  largo. 
Computed  profiles  for  other  values  of  F  and  0  are  qualit  itively  similar  to  these  two.  We  were  able 
to  find  solutions  for  all  the  values  of  F  >  1  which  we  tried. 

For  F  <  1,  (2.6)  has  purely  imaginary  roots  and  then  (2.5)  corresponds  to  a  train  of  small 
amplitude  waves  at  infinity.  Although  we  did  not  compute  them,  we  expect  that  there  are  solutions 
for  F  <  1  with  such  wavetrains  at  infinity. 


S.  Pouring  flows  and  weir  flows 

We  now  consider  flows  with  two  free  surfaces  over  a  thin  wall  in  water  of  infinite  depth.  (See  Figure 

1. )  The  flow  in  Figure  1  with  arbitrary  wall  slope  0  is  a  generalization  of  the  flow  with  0  =  x/2 
calculated  by  Vanden-Broeck  and  Keller  [l].  Solutions  for  arbitrary  0  can  be  obtained  by  using  the 
procedure  described  in  Section  3  of  their  paper  with  their  equation  (3.8)  replaced  by 

C  =  -e<(f-/,)(l  +  t)J*^[-lnc(l  +  tl)]1/3exp(f;j;ntn)  (5.1) 

nsO 

For  0  =  r/2,  (5.1)  reduces  to  their  equation  (3.8).  A  computed  profile  obtained  by  using  (5.1) 
with  0  =  x/3  is  shown  in  Figure  1.  Similar  results  obtained  with  different  values  of  0  indicate 
that  there  is  a  unique  flow  with  two  free  surfaces  for  each  value  of  0  in  0  <  0  <  x/2.  However  the 
analysis  presented  in  the  previous  sections  suggests  that  there  are  in  addition  solutions  with  one 
free  surface. 

In  order  to  compute  these  extra  solutions  we  consider  the  flow  configuration  shown  in  Figure 

2.  We  choose  0  =  0  on  the  wall  ICJ  and  <p  =  0  at  the  point  C.  Let  Q  be  the  value  of  rfr  on  the 
free  surface.  We  introduce  dimensionless  variables  such  that  Q  =  g  =  1.  In  these  variables,  the 
Bernoulli  equation  yields 

(Vy>)2  -  2*  =  0  on  tfr  =  1.  (5.2) 

The  complex  potential  plane  is  shown  in  Figure  5a.  By  using  the  transformation  (2.10)  as  before, 
we  map  the  flow  domain  into  the  interior  of  the  unit  circle  in  the  t-plane  shown  in  Figure  5b.  It 
can  be  shown  that  the  following  expression  is  a  suitable  expansion  for  the  complex  velocity: 

C  ~  e<<t-«(l  +  f)^r,[-lnc(l  +  t)],/3exp  (  £  UntA  (5.3) 

'bsO  ' 

Here  c  is  an  arbitrary  constant  between  0  and  0.5.  As  before  we  choose  e  —  0.2.  The  coefficients  Un 
are  determined  by  the  collocation  procedure  of  Section  2  with  the  parameter  F  in  (2.16)  replaced 
by  1.  A  typical  profile  for  0  =  x/3  is  shown  in  Figure  2. 
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